GET P, N 



= 0 



i 



|0 



1 = 1 + 1 

G = GCD ( P, THOUSAND.PRIMES [I]); 



<^T>U 



GETP, M,E 



DETERMINE M* 1 MOD P 
BY CALCULATING ' 
GCD (P, M) 



RETURN FALSE 



i 



APPLY SIGNED DIGIT 
ALGORITHM TO E 



i 



|6 



1 i 







RETURN TRUE 





CALCULATE M E MOD P 



FiC 3 



Pi4 i 



GETU.V 



U3 = U; V3 = V; 
U2 = 0;V2 = 1; 



i 



Zio 



V3 = 0 



T 





RETURN 




U3, U2 



Q = LU3/V3j; 
T3 = U3-Q*V3 
T2 = U2-Q*V2 
U3 = V3; V3 = T3 
U2 = V2; V2 = T2 



GETU.V&M 



U3 = U;V3 = V; 



AY 1 



D = 8 4-<M = 32>— ► 
-p 1 y > ' ri 



D = 5 



4ft 



[A,B] = TOP2M(U3 I V3, M); 




▼ S 



L0G(U3) - L0G(V3) > D 




Q1 = LA/Bj; 
T3 = U3-Q1 *V3 



m4 



K 




CALCULATE 
XA + YB<2 M 




> A 


T3 = XU3 + YV3; 
U3 = V3;V3 = |T3| 



[A, B] = TOP2M(|T3|, V3, M); 
Q2 = Q1- (LA/Bj+1); 



■A jr3l>V3 > 



t4 



4£7 



Q2 = Q1-1 



T3 = U3 - Q2 * V3 



N 



43P 



Q = l_U3/V3j 
T3 = U3-Q*V3 



U3 = V3;V3 = T3; 



— 



43* 



Y = 1 



Y = GCD (U3, U MOD U3); 
Y = GCD (Y, V MOD Y); 



RETURN Y 



GET A, B, M < 




r 



^1 



6\y 




M = 32 




■6\o 



Y 



[X, Y] = MOD_EUCLID_GCD(A, B) 



[X, Y] = MOD_LEHMER_GCD(A, B) 



* X. 



RETURN X,Y 



FIG. 5 



GET U3, V3, M 



PB1 = LOG (U3) 



BP2 = BP1 - 2M + 1 



RETURN U3[BP1:BP2], V3[BP1:BP2] 



FIC4A 



GET U, V 



U3 = U;V3 = V; 
U2 = 0;V2 = 1; 



U>lO 



Y 




L0G(U3) - LOG (V3) > 32 




[SU,SV] = TOP2M(U3 > V3 l 16); 
D = LOG (U3)-32; 
J = 0; DONE = FALSE; 
X[0] = 1;X[1] = 0;Y[0] = 0,Y[1] = 1; 



LP2.4 



Q = LSU/SVj; 
X[J+2] = X[J]-X[J+1]; 
Y[J+2] = Y[J]-Y[J+1]; 
T = SU-Q*SV;SU = SV;SV = T; 
J=J+1; 



DONE = SV<-Y[J+1]OR 
SU-SV<X[J+1]-X[J] OR 
SU<2 D ; 



J EVEN 



bl-2- 



DONE = SU<-X[J+1]OR 
SU-SV<Y[J+1]-Y[J]OR 
SU<2 D ; 



T 



DONE = TRUE?> 



T3 = X[J-1] * U3 + Y[J-1] * V3; T2 = X[J-1] * U2 + Yfl-1] * V2; 

V3 = X[J] * U3 + Y[J] * V3; V2 = X[J] * U2 + Y[J] * V2; 
U3 = T3;U2 = T2; 



Q = LU3/V3j; 
T3 = U3-Q*V3;T2 = U2-Q*V2 
U3 = V3; V3 = T3; U2 = V2; V2 = T2 



F1&. U> 



U1 = (U3-U2*V)/U; 
RETURN U1.U2 



GET U, V 



7 I O 



i 



U3 = U;V3 = V; 
U2 = 0;V2 = 1; 



< ^U3<2 16 ^ > 

i 



7'*/ 



7/6 



U1=(U3-U2*V)/U; 
RETURN U1.U2 



Q 


= LU3/V3j; 


T3 = 


U3-Q*V3; 


T2 = 


U2-Q*V2; 


U3 = 


V3;V3 = T3; 


U2 = 


V2;V2 = T2; 



-7(6 



GETU.V&M 



WO 



U3 = U;V3 = V; ^i-p 
U2 = 0;V2 = 1; rZ> 



131 



(A, B] = TOP2M(U3, V3, M); 




L0G(U3) - L0G(V3) > D 



Y 






Q1 = LA/Bj; 


1 — ► 


T3 = U3-Q1 *V3 











CALCULATE 
XA + YB < 2 M 



1 



y 



?3b 



Q = Q1 



ri 



>V3 




[A, B] = TOP2M(|T3|,V3, M); 
Q2 = Q1- (LA/Bj + 1) 



Q2 = Q1-1 



T3 = U3-Q2*V3 



T3 = XU3 + YV3; 
T2 = XU2 + YV2; 



T2 = -T2 



Q = Q2 



4-<T37o> 



84! 



^KT3<0>^ 



U3 = V3;V3 = |T3| 
U2 = V2;V2 = T2 



Q = LU3/V3j 

T3 = U3-Q*V3 




T2 


= U2-Q*V2; 




113 


= V3;V3 = T3; 




U2 


= V2;V2 = T2; 





- c± < v V3 = 0^> 



Y=1 



Y = GCD (U3, U MOD U3); X = LU3/Y.I; 

COMPUTE X" 1 WRT U BY GCD (U, X)- 
U2 = (U2*X" 1 )MODU 
Y = GCD (Y, V MOD Y); 



RETURN Y, U2 



